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The dynamical properties of road traffic time series from North- Rhine Westphalian motorways are 
investigated. The article shows that road traffic dynamics is well described as a persistent stochas- 
tic process with two fixed points representing the freeflow (non-congested) and the congested state 
regime. These traffic states have different statistical properties, with respect to waiting time distribu- 
tion, velocity distribution and autocorrelation. Logdifferences of velocity records reveal non-normal, 
obviously leptocurtic distribution. Further, linear and nonlinear phase-plane based analysis methods 
yield no evidence for any determinism or deterministic chaos to be involved in traffic dynamics on 
shorter than diurnal time scales. Several Hurst-exponent estimators indicate long-range dependence 
for the free flow state. 

Finally, our results are not in accordance to the typical heuristic flngerprints of self-organized criti- 
cality. We suggest the more simplistic assumption of a non-critical phase transition between freeflow 
and congested traffic. 

PACS numbers: PACS. 45.70.-n Granular systems; traffic flow. 



I. INTRODUCTION 

Traffic flow prediction, particularly in close connection 
with the avoidance of jams, is a challenging, yet hitherto 
unreached target. 

Until now, due to restricted access to records, simu- 
lation models provided the predominant approach to 
understand traffic dynamics. Several approaches have 
been developed which are based on partial differential 
equations (0,0) j or cellular automata models as the 
widespread Nagel-Schreckenberg model (^]). A com- 
prehensive overview of results from time series analysis 
from real traffic records was published by Q. In earlier 
research on the database that our study relies on, diur- 
nal, weekly and annual cycles in traffic density as well as 
velocity was reported in details by 0. Autocorrelation 
and time-headways of traffic records are demonstrated 
to vary state-dependently (0,0), congested traffic 
revealing a more persistant autocorrelation. 
Intuitively, traffic dynamics conforms rather to a 
stochastic than deterministic(-chaotic) process. A 
rigorous statistical inference however, to the best of our 
knowledge has not yet been achieved. 

This paper is organized as follows: We first intro- 
duce the dynamical phase-plane reconstruction from traf- 
fic records by fundamental diagram and delay-plot, to 
point up that traffic dynamics consist of two heteroge- 
neous states. The further analysis focuses on separated 
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sections of either free flow or congested traffic regimes. 
We then turn to phase-plane based methods such as cor- 
relation integrals, and surrogate based local linear pre- 
dictions to demonstrate that traffic dynamics on below 
diurnal time scales has a predominantly stochastic na- 
ture. 

Long-range dependence is tested from several measures. 
To exclude possible effects of nonstationarity, the latter 
measure is compared with appropriate phase randomized 
surrogates. Nonlinearity will be discussed by application 
of the surrogate based time-reversibility test. 



II. DATA ANALYSIS 
A. Methods 

1. Phase-randomized surrogates 

In time series analysis, phase-randomized surrogate 
(PRS) time series (Q) can be applied as a version of 
bootstrapping to clarify and quantify statements about 
the presence of nonlinear effects. PRS series reveal the 
same linear statistical properties as their original and can 
be produced at will. Possible nonlinearities, as nonlinear 
determinism beyond the autocorrelation of the original 
time series will not be reproduced by their surrogatiza- 
tion, or changed by interpretation as a spectral property. 
In summary, PRS time series are produced by multiply- 
ing the Fourier-spectrum of the original records with ran- 
dom phases and hereafter performing a backtransforma- 
tion (for details see or [loj). 
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2. Nonlinear methods 

In this paper we will make use of linear and nonlinear 
phase-plane based measures such as correlation dimen- 
sion and local linear prediction. Such methods are usu- 
ally applied to time series with the intention of identify- 
ing the presence of nonlinear, possibly chaotic dynamics. 
Since it is hardly possible to formally prove the absence 
of any deterministic property, we intend to point out this 
absence by comparing (nonlinear) statistics for original 
data vs. their appropriate surrogate substitutes. 

B. Records 

Freeway traffic in North-Rhine Westphalia (Germany) 
is continuously monitored at approximately 1400 road 
locations by means of built-in loop detectors. For every 
appearance of a vehicle these detectors record: 

1. time, 

2. velocity, 

3. type of vehicle, 

4. length of the vehicle. 

This study is based on two different types of static loop- 
detector recordings: 

1. Single car records: 

Only a few exceptional time series have been 
recorded with a notebook PC attached to the loop 
detectors, 

2. minute-aggregated data: 

These data are obtained from the same loop detec- 
tors as single- car data. However, instead of im- 
mediate recording, the samplings are exponentially 
smoothed and aggregated in 1-minute intervals. 

Both single-car and minute aggregated records are 
coarse-grained, since all records are denoted in [ "0" ... 
"254" ], while "255" denotes faulty results. Due to their 
higher resolution, single-car data provide rare, but the 
most detained (and unspoilt) information, particularly 
for short time scales. 

For more details of sampling and processing read Q and 
0. Both articles provide a detailled introduction into 
practial aspects of road traffic data. 

III. RESULTS 

A. Time course of traffic dynamics 

Fig. n presents a section of typical single car freeway 
traffic records. In a), the velocity time series appears 
more or less regularly fluctuating, except for occasional 
abrupt drops in velocity, (b) is a one-day sequence of 
single-car records comprising a jam episode. 
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FIG. 1: velocities, 

a) minute aggregated records, covering 60.000 min ~ 1000 h 
^ 40 d. 

(b) 1 day of single car data comprising a jam event, arrows 
indicate sections of jam and freeflow traffic state that will be 
analyzed in the following. 
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FIG. 2: Fundamental diagram: plot of occupancy [sum of 
vehicle lengthes per road kilometer] vs. flux [vehicles per 
minute] for minute aggregated traffic records. 



1. Fundamental diagram 

In the context of traffic analysis, the fundamental dia- 
gram, well known to applied scientists, denotes the plot of 
flux vs. occupancy, in most cases graphed from smoothed 
model output data (e.g. [3) . 

Fig. 121 was graphed from minute aggregated records. Due 
to the discreteness of the latter, in a plot like Fig. |21 
some hundreds of thousands of data would fall into a few 
bins. To improve the visualization, we added uniformly 
distributed independent random noise [—0.5 < ^ < 0.5] 
(the noise level scales below the resolution of the signal) 
to the data. Fig. □ gives an impression of traffic dynam- 
ics, that undergoes transitions between two attractive re- 
gions representing freeflow and jammed state. Whereas 
the freeflow regime (high velocity, low occupancy) is sit- 
uated transversally on the left hand side, the congested 
state associates with a larger realm of points in the center 
and on the right hand side. 
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FIG. 3: 

(a) Delay-plot minute aggregated velocity data covering one 
month, r = 1 minute. 

(b) Delay plot of single-car data - (1 day), r = 25 sec. 



1. Local stationarity assumption 

Naturally, the double fixed point structure, visualized 
in figures |21 and |31 gives a strong indication against sta- 
tionarity for the overall process of traffic dynamics that 
comprises two traffic states. In the following we will 
therefore constrain the analysis to selected sections of 
either freeflow or jammed traffic that are indicated by ar- 
rows in Fig. ^ To apply methods that require regularly 
sampled data we transform these sections of single-car 
data equidistant by aggregation and linear interpolation, 
expecting that this procedure does not have substantial 
influence on the results. 



C. Distribution of intervals between consecutive 
events (time-headways) 

Time-headway distributions from single car data have 
already been reported in for different traffic states. 
According to our results they reveal an approximately 
lognormal distribution with different parameters in de- 
pendence of the traffic state (Fig0] (a) and (b)). For 
freeflow traffic the Kolmogorov-Smirnov test statistics 



B. Delay plot 

Fig-Elcan be interpreted as "phase plane" formed from 
records. An alternative, more practicable method to ob- 
tain a comparable clue on phase plane is delay coordi- 
nate embedding, which denotes a n-dimensional plot Xt 
vs. Xt+r vs. ... vs. ^ of a time series 

Xt,t = l,2,...N. 

The well-known general results by Takens ^3 state that 
the dynamics of a system recovered by delay coordinate 
embedding are comparable to the dynamics of the origi- 
nal system. A low dimensional deterministic-chaotic at- 
tractor thus can be graphed from each of its observed 
variables as a topologically equivalent structure to what 
one would obtain from the graph of its variables in a 
sufficiently dimensioned delay plot. Since there is no 
straightforward way to determine which dimension is suf- 
ficiently large, several dimensions need to be examined. 
According to [lOj an optimal delay r approximately cor- 
responds to the empirical autocorrelation function (ACF) 
at r(r) — 1/e. Fig.|21i) shows a delay plot of the veloc- 
ity series Xt used in Fig. |21 here plotted in time-delayed 
coordinates x(t — t), x{t) and x{t + t), t > denoting 
the delay-time. In Fig. O (a) and (b) the data scatter 
around two condensed regions, which can be identified 
as congested and freefiow traffic. 

We present Fig. |3| to visualize a clearer two fixed-point 
structure than in the "traditional" plot Fig.|5| Moreover, 
though the single car data base is not sufficient to obtain 
a fundamental diagram. Fig. |21 (b) gives an indication of 
comparable dynamics in single car data. 



D = = 0.0000804 (1) 

n ~ 500 

performs below the tabelled value 

1 949 

i?„^o.ooi - , • = 0.087. (2) 
V?^ = 500 

Thus, this test on distributional adaptation does not 
state the rejection of the null-hypothesis of lognormal dis- 
tribution. In Fig. 01a) however, a deviation in the right 
wing (reminding to a fat tail) is observed. The finite 
left tail of the distributions probably reflects the neces- 
sity to keep a security distance between vehicles. Data 
of jammed traffic (Figure ^b)) are comparably scarce. 
Little, if anything, can be inferred from them. 



D. Self Organized Criticality 

Previous authors (fT5|) already suspected that road 
traffic has a selfsimilar nature in the context of the Self- 
Organized Criticality (SOC) models. According to such 
models, increasing traffic load would produce a "critical" 
situation, that, at its critical point, occasionally relaxes 
catastrophically (e.g. as sandslides in the sandpile model 
[T^). Close to the critical point, such a system generates 
power law behaviour, observable in leptocurtic distribu- 
tions, slowly converging variance, lack of characteristical 
scales and 1/f noise. 
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FIG. 4: (a) Time interval histogram of entire single car data 
series (dashed line), fitted lognormal distribution (solid line)), 
(b) Time interval histogram during jammed state, fitted log- 
normal distribution (dashed line). 
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FIG. 5: Histogram of single car velocities of 12 different 
single-car highway traffic data sources (solid line), 
interpretation as addition of 2 Gaussian distribution curves 
(dashed line). 




FIG. 6: (a) Histogram of differenced single car velocities se- 
lected for jam and freeflow epsidodes, 

(b) cumulated distribution of logdifferenced single car veloci- 
ties selected for left(l) and right (r) wings. 



1. Distribution of velocities and velocity differences 

Fig. [SI shows the histogram of traffic velocities of all 
single-car data (12 different locations). In some loca- 
tions slow congested traffic and jams appear as a smaller 
second peak that, e.g. for smaller data quantity, could 
be misinterpreted as fat tail in the low-speed end of the 
empirical velocity probability distribution function. 
Comparable to the well-known heavy-tailed distributions 
of logdifferenced financial time series, in fig. (a) we 
observe a clearly non-normal distribution in differenced 
velocity records, separated for either jam and free-flow 
records. This holds for logdifferenced data as well (not 
shown here). The histogram looks more leptocurtic for 
free-flow than for jammed traffic records. 
The plot of the cumulated distribution function in 
double-log coordinates provides a clue if the asymptotic 
behaviour of the functional form of the cumulative dis- 
tribution is "visually" consistent with a power law, 

P{G >^*)^^ (3) 

where a is the exponent caracterizing the power law de- 
cay, 

G{t,M)=\n{xt+^t)-Hxt) (4) 

Together with other indications, such a power law 
can be regarded as a feature which is characteristic for 
SOC processes ([3). 

Conversely, the cumulated distribution of differenced ve- 
locities, separated for right and left wings of either jam 
and freeflow records, displays no clear scaling region. 
Non-normal distribution as well as lack of scaling is also 
observable for the larger database of minute aggregated 
records (not shown here). 

E. Long-range dependence 

Scientists in diverse fields observed empirically that 
positive correlations between observations which are far 
apart in time decay much slower than would be expected 
from classical stochastic models. In time scries such cor- 
relations are characterized by the Hurst exponent H. 
They are often referred to as Hurst effect or long-range 
dependence (LRD). 0.5 < H <1 reflects long-range posi- 
tive correlations between sequential data. H = 0.5 corre- 
sponds to sequential uncorrelatedness (known for white 
noise). Brownian motion, the trail of white noise, is char- 
acterized hy H = \. 

Since long-range dependence (LRD) is defined by the au- 
tocorrelation function (ACF), theoretically, the shape of 
the ACF provides an indication for LRD in road traf- 
fic. For LRD series, the ACF at large lags should have a 
hyperbolical shape: 

r(r) cx r^^^^T ^ cx) (5) 
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(na)- 

The practical ability to assure an algebraic decay of the 
ACF however is low, making such an approach inviable 
for data analysis. For comparable reasons, from the tail 
of the distribution, additional information is hardly ob- 
tainable; statistics here are generally poor (18]). The 
discreteness of car traffic data additionally diminishes the 
quality of such estimations. 



1. Hurst-exponent estimation 

The estimation of the Hurst-exponent (H) from em- 
pirical data is not a simple task. Several studies (e.g. 
|l7j.p^.[20l|) estimate the Hurst exponent H from differ- 
ent measures. 

Synthetically generated fractional Brownian motion or 
fractional ARIMA (autoregressive integrated moving av- 
erage) series are characterized by a generalized (or global) 
H . Such so called monofractional series are known to re- 
veal fluctuations on all time scales. 

They will produce unambiguous evidence for fractional- 
ity, whereas a more general class of heterogeneous signals 
exist that are made up of many interwoven sets with dif- 
ferent local Hurst exponents, called multifractional f[2ll|). 
It is a frequent experience, that graphical methods to 
test for LRD show no clear scaling for such series. Our 
own experience is, that weighted sums of synthetically 
generated random walks with different characteristical 
scales may as well reveal straight fractional scaling in 
some plots, as crossover behaviour according to other 
methods. Furthermore, some methods of i/-estimation 
sensitively depend on the distribution of the data. 
The main criticism against i7-estimates is based on the 
experience that instationary data may, at least in some 
cases, produce estimates that erroneously indicate frac- 
tionality. Thus, we are interested in the robustness of 
i/-estimators, if possible effects of instationarity are ex- 
cluded. Phase-randomized surrogates (PRS) based on 
original traffic records are random sequences with the 
same first and second order properties (the mean, the 
variance and the auto-covariance data, but which are oth- 
erwise random. 

Since fractionality is a spectral property, and PRS fully 
recover the latter, i?-estimation of PRS hence provides 
an approach to exclude possibly misleading effecs of in- 
stationarity, albeit not to differentiate between monofrac- 
tional and heterogeneous signals. 

To obtain reliable H-estimates despite possible effects 
of nonstationarity, we apply a variety of the most famil- 
iar methods to jam- and freeflow traffic records. For a 
detained discussion about the application of methods to 
conclude on LRD, e.g. aggregated variance method, R/S 
plot, periodogram method and wavelet-based Whittle- 
estimator on nonstationary data read p^ . Detrended 
fluctuation analysis (DFA) (^2) denotes the root mean 




FIG. 7: 

R/S pox plot of (a) freeflow and 

(b) congested velocity records, approximation of scaling by 
Hurst exponents (a) H = 0.71, (b) H — 0.81 (solid lines). 
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FIG. 8: 

(a) Spectrum from selected freeflow series (dots), exponen- 
tially weighted moving average (solid line), scale exponents 
a « 0.71, (dashed line). 

(b) Spectrum from selected jam series (dots), exponentially 
weighted moving average (solid line), scale exponent a = 0.88 
(dashed line) 



square fluctuations 
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around least squares line trend fits yn{k) for equal box 
length n of the integrated series 



yik) = ^{xt - x), 



N 



(7) 



A straight line in the double logaritmic plot that indi- 
cates scaling has the slope 2H . The method should pro- 
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0.611 




0.875 




DFA 


0.685 




1.107 




" for PRS 


0.669 


0.005 


1.232 


0.017 



TABLE I: Hurst-exponent estimators from traffic records for 
different methods: R/S : rescaled range analysis, 
a.v.: aggregated variance method, 
a.a.: aggregated absolutes, 
spc: graphical estimation from the spectrum, 
DFA: detrended fluctuation analysis. PRS denotes the appli- 
cation of the above methode to phase-randomized surrogates, 
(T denotes the standard deviation. 



vide robust estimates even for nonstationary time series. 
Table ^ displays the results. 

We also applied the wavelet-based Whittle estimator 
(01). Despite its postulated robustness against insta- 
tionarity, and despite iJ-estimates that compare to ta- 
blcQI we do not show the graphs here, since, particularly 
for the jam series, the wavelet-spectrum offers to many 
possibilities of parametrization, as. e.g. the choice of the 
wavelet function, octaves etc.. 

Fig. [3 presents the R/S pox plots of freeflow and jam 
records. Particularly for freeflow data exact straight scal- 
ing is not observable. The same, even more, holds for 
Fig. |Sla). In anao gy t o the the modified periodogram 
method outHned in [13, the logaritmically spaced spec- 
trum was divided into 60 boxes of equal length. The least 
squares ft was performed to averages of the data inside 
these boxes, to compensate for the fact that most of the 
frequencies in the spectrum fall on the far right, whereas 
for LRD-investigation the low frequencies are of inter- 
est. Fig. IHIa) gives the strongest indication that traffic 
dynamics can not be characterized as monofractional as 
most of the common Hurst-cstimators would indicate. 



F. Time reversibility test 

An important property to differentiate between linear 
and nonlinear stochastic processes is time-reversibility, 
i.e. the statistical properties are the same forward and 
backward in time. From this test one can not judge 
whether the data correspond to any ARMA-model, since 
theoretically, time asymmetry might be caused by non- 
Gaussian innovations. Apart from on-ramps, traffic dy- 
namics on short time scales anyway is unlikely to be sub- 
stantially influenced by external noise. 
The following expression is outlined as a measure to con- 



clude on time reversibility of time series (| 



Qij) 



Xt) 



E[{xt+r - Xt)^] 



(8) 



wherein r denotes the delay time and E represents the 
time average. The basic idea behind it is to compare the 
time reversibility test statistics of the original data Q{t) 
with confidence bounds from corresponding test statistics 
QsuriT), generated from surrogate series: 



^(C^,r) < Q{t) < Qsurr{Ca,T); 



(9) 



for some critical ; Cq . 

The results for a surrogate-based test are usually re- 
ported as significances: 



Sir) 



ViQjr) - {Q{r)),urry 

surr 



(10) 



where: 

Q{t) test statistics, 

{Q{T))surr mean, 
<j{Q{T))surr Standard deviation. 

The test is based on the assumption that the surrogate 
test statistics for a given lag are approximately Gaussian 
distributed. 

The statistical properties of the examined surrogate time 
series are the same forward and backward in time. Thus 
they comply with the null-hypothesis of time-reversibility 
which will be tried to reject by the test. 
The evaluation of significances for more than one lag 
leads to the statistical problem of multiple testing. This 
has severe implications on the probability to reject the 
null- hypothesis. The Bonferroni-correction of the signif- 
icance level must be taken into account: 



d = 1 - (1 - a)" 



(11) 



wherein n denotes the number of independent tests. 
Practically, Bonferroni- corrected confidence bands ren- 
der little diagnostic power to detect a violation of the 
null- hypothesis. A corrected significance level, for ex- 
ample, 1 — d = 0.95 for 100 independent tests requires 
1 — a « 0.9995. In most cases, however, Q{t) is au- 
tocorrelated to an unknown extent, what diminishes the 
number of independent tests and, for rejection of the null- 
hypothesis, results in a conservative test design. 
In Fig. El surrogate-based time-reversibility tests for jam 
and freeflow traffic states are graphed. Under the as- 
sumption of 100 independent tests for a = 0.001, the 
corrected significance level is: d = 1 — (0.9999)"^^°" w 
0.905, which, though not acceptable as safe statistical 
inference, gives a vague information that freeflow traf- 
fic dynamics is more likely time-irreversible than time- 
reversible. For the jammed state the observed 20 devi- 
ations of the confidence bands gives a comfortably safe 
rejection of the nuU-hyothesis, particularly for short time 
scales, but also for larger r. Since, even for the naked eye. 
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FIG. 9; (a) Q(r) from freeflow series (crosses), lines indicate 
confidence limits from 500 surrogate realizations, (b) Q{t) 
from jam traffic series (crosses), lines indicate confidence lim- 
its from 50 surrogate realizations. 

the test statistics is substantially correlated, the test pro- 
vides a safe rejection. For the freeflow state, 14 deviations 
from Hq axe also statistically indicative, albeit less corre- 
lation among the test statistics is observable. Both traffic 
states thus are likely to reveal time-irreversible statistical 
properties. 



G. Recurrence plot 

For a time series Xt the recurrence plot is a two- 
dimensional graph that is formed from embedded vectors 

y{t) = {xt,Xt+r, ■ • ■ , Xt+(^E-l)T) (12) 

for embedding dimension E and lag r. These vectors are 
compared if they are in e-proximity of another y{t + At). 
If 

\\y{t)-y{t + m\<e, (13) 

a black point is drawn at {t, At). For each e, r, m (with: 
m the embedding dimension, r the time lag, e the vari- 
able error distance) an individual recurrence plot is ob- 
tainable. 

Since the differences 

yiU) - y{tj) = = Cj- , = y{tj) - y{ti) (14) 

are identical, the plot consists of two symmetric triangu- 
lar graphs along a black (since i = j) diagonal line. 
Except for horizontal and vertical stripes (that might re- 
flect temporal (auto-) correlations), the recurrence plot 
of freeflow traffic Fig. 1101 is very much in remedy of what 
one would observe for a recurrence plot of a white noise 
series. 




FIG. 10: Recurrence plot of freeflow data, r = 1.7 seconds, 
e = 0.01, TO = 3. 

H. Correlation integral 

There exist 

Nt = ^N{N - 1) 

independent radii c„ (since c„_m — Cm,n)- The density 
of a recurrence plot as functional of e 

t=N n-1 

C{e, D, t) = 1/Nt 5m ^(e - - ^t-r\E) (15) 

t = l T=l 

{6 denoting the Heavyside step function, 

with e{z) = 1, for z > 0, e{z) = for z < 0) 

is called the correlation integral. The resulting Cr{e) is 

sketched in a double logarithmical Grassberger-Procaccia 

plot ([l^l) in dependence of e. 

The correlation integral is plotted for varying dimension 
as well as varying e. If a noise-contaminated determin- 
istic process is regarded, from a sufficient embedding di- 
mension, parallel slopes for varying dimensions indicate 
power law scaling in a region which is situated above a 
certain e that represents the noise range. 

Figure shows Grassberger-Procaccia plots for (a) 
free flow traffic and (b) jammed state records for error 
distances e of values 0.1 . . . 100 and embedding dimension 
2 ... 25. Both graphs fail to reveal any scaling region, 
moreover there is obviously no difference to Grassberger- 
Procaccia plots of appropriate surrogate realizations. 
The dimensions of merely stochastic systems appear in- 
finite, therefore for this case it is a typical result, that 
the correlation integral reveals the embedding dimension 

(H). 

1. Casdagli test 

The local linear prediction of a time series in delay rep- 
resentation Xt is achieved by determination of a matrix 
A that minimizes the prediction error: 

't'= E {3Lt+i-Ax^^h)\ (16) 
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FIG. 11: 

a) Grassberger-Procaccia plot for (a) freeflow- and (b) 
jammed state traffic data for varying dimensions and lag t = 7 
seconds (solid lines). Grassberger-Procaccia with identical 
parameters for appropriate surrogate series realizations are 
plotted in dashed lines. 

where Ut denotes the e— neighbourhood of xt excluding 
Xt itself. In some analogy to linear regression the predic- 
tion is: 

x;+i=Axt + bt. (17) 

Local linear models are suggested a test for nonlinearity 
( 26] ). The average forecast error is computed as func- 
tion of the neighbourhood size on which the fit is per- 
formed. If the optimum occurs at large neighbourhood 
sizes, the data are (in this embedding) best described by 
a linear stochastic process, whereas an optimum at rather 
small neighbourhood sizes supports the idea of existence 
of a nonlinear, almost deterministic, equation of motion 

(H). 

In Fig. El the Casdagli-test was performed for origi- 
nal as well as for surrogate data. In the context of the 
Casdagli-test it is meaningless some of jam records fall 
not within the surrogate-based confidence intervals. The 
qualitative comparability gives an indication for predom- 
inant stochasticity particularly in freeflow traffic records. 

IV. DISCUSSION 

According to our analysis of traffic records, traffic dy- 
namics is a two fixed point stochastic process, while the 
fixed points reflect the jam and freeflow regime. 
The abrupt transitions between the traffic states imply 
nonlinearity in the overall traffic dynamics. A variety 
of methods, more or less sensitive towards nonstationar- 
ity, yields Hurst-exponents that indicate long-range de- 
pendent dynamics in particular for freeflow traffic. Dif- 
ferenced as well as logdifferenced velocity records reveal 
heavy tailed distribution, however for both there is no 
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FIG. 12: Casdagli plot of (a) free-ffow and (b) jammed state 
traffic data (crosses), compared to confidence bands, gener- 
ated from to PRS series (dotted lines). 



clear scaling region observabele to estimate scaling expo- 
nents. 

From our results it must be concluded that below diur- 
nal time scales traffic data in jammed as in freeflow state 
exhibit neither deterministic nor low dimensional chaotic 
properties. 

The main intention of this article is to outline an overview 
of the stochastical properties achieved by data analysis 
of single-car road traffic records. 

Attending the problem of criticality in road traffic 
records, we find that as well the two fixed-point dynam- 
ics as the distribution of (differenced) velocities are con- 
trary to the typical features of processes governed by self- 
organized criticality. This lets us rather suspect the rise 
of jams in the context of a (eventually, but not neces- 
sarily, critical) phenomenon linked to a phase-transition. 
For such a model hypothesis, known e.g. in equilibrium 
thermodynamics, the point of transition can be reached 
by fine tuning of a parameter. This must be distinguished 
from self-organized criticality, which represents the clas- 
sification for systems attracted permanently by variable 
critical states. 

Contrary to the well-known conceptual analogy between 
traffic and granular flow, we rather propose an intu- 
itive analogy of traffic dynamics with the condensation of 
steam to water. In contrast to the condensation of water 
driven by withdrawal of heat, free flow traffic condenses 
to higher particle density by an increase of trafficants in 
this picture. This increase can be interpreted as rised 
pressure. In accordance to such considerations and the 
empirical results of 4] increasing traffic load (or input to 
the motorway) produces a (in the more popular sense) 
"critical" tension that relaxes in an abrupt transition to 
a jam. In this instructive example, traffic accidents, con- 
struction sites or slow vehicles could act comparable to 
condensation cores by exerting strong nonlinear negative 
feedback on the upstream traffic. The fine tuning pa- 
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rameter thus is the capacitiy of the motorway, limited by 
traffic load, accidents or construction sites. 
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